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Abstract 

The multi-component nonlinear Schrodinger equation related to C.I ~ Sp(2p) /U(p) and 
D.III ~ 50(2p)/(7(p)-type symmetric spaces with non-vanishing boundary conditions is solv- 
able with the inverse scattering method (ISM). As Lax operator L we use the generalized 
Zakharov-Shabat operator. We show that the ISM for the Lax operator L(x, A) is a nonlinear 
analog of the Fourier-transform method. As appropriate generalizations of the usual Fourier- 
exponential functions we use the so-called "squared solutions", which are constructed in terms 
of the fundamental analytic solutions (FAS) x (a;, A) °f L(x,X) and the Cartan-Weyl basis 
of the Lie algebra, relevant to the symmetric space. We derive the completeness relation for 
the "squared solutions" which turns out to provide spectral decomposition of the recursion 
(generating) operators A± , a natural generalizations of \ -4- in the case of nonlinear evolution 
equations (NLEE). 

1 Introduction 

The integrability of the scalar nonlinear Schrodinger equation (NLS) with vanishing boundary 
conditions (v.b.c): 

iqt + qxx + ^qix^qix^t) = (1.1) 

was discovered by Zakharov and Shabat in their pioneer work |21j. Soon after [22 Zakharov and 
Shabat proved the integrability and the physical importance of the NLS with constant boundary 
conditions (c.b.c): 

iq t + 2q xx -2(\q(x,t)\ 2 - p 2 )q(x,t) = 0, lim q(x,t) = q±, (1.2) 

x — >±oo 

where the asymptotic values q± satisfy |<7±| 2 = p 2 . Notice the sign difference in the cubic nonlin- 
earity as well as the additional term with the chemical potential p. 

Both versions of NLS equation served as models on which generalizations were made. The 
simplest non-trivial multicomponent generalization of NLS is the vector NLS known as the Manakov 
model [IS] : 

i~qt + ~q X x + 2(~q^~q(x, t))~q(x, t) = 0, (1.3) 

where ~q(x,i) is an n-component complex-valued vector vanishing fast enough for x — > ±oo. The 
c.b.c. version of vector NLS 

i~qt + ~q xx ~2 ((~q^~q(x, t)) - p 2 ) ~q{x, t) = 0, (1.4) 



where lim x ^± OQ ~q(x, t) = ~q± and q _ = Uo~q+ where Uo is constant unitary matrix also finds 
applications. Here p 2 = ~~q±~q±. 

Equations ( |1.1[ ) and ( 1.3 1 are particular cases of matrix NLS which is obtained from the system: 



2qrq(a;, t) = , -ir t + r xx + 2rqr(x, t) = 0, 



(1.5) 



after imposing appropriate involution (reduction) compatible with the evolution of (1.51. Here q 
and r are n x m matrix-valued functions of x and t. One such involution is: 



q' u _\ B± = diag(ef,...,e±), (ef ) z = 1, 
and the corresponding MNLS acquires the form: 

i q t + q ra + 2qB_q t 57^ = 0, 



(1.6) 



(1.7) 



For n — rn — 1 and r = q* the system goes into the scalar NLS ( 1.1 ); for m — 1 and n > 1 and with 
appropriate choice of involution ( 1.6 1 the system is transformed into the Manakov model ( 1.3 1. All 
these versions are solvable with the ISM. The ISM is applicable to nonlinear evolution equations 
(NLEE) if they can be represented as compatibility condition of two linear problems [TH1 123 12 H] : 



[Z(A),M(A)]=0, 



(1.8) 



which holds identically with respect to the spectral parameter A. 

The two linear operators L(X) and M(X) in the Zakharov-Shabat system (Z-Sh) for the MNLS 
on symmetric spaces associated with the simple Lie algebra g ~ C r and g ~ D r with (v.b.c.) are: 



Lip = 
Mil) = 

Q(x,t) = 



d 

i— + Q(x, t) - A <t 3 ] ip(x, t, A) = 0, 



ig- t +V 2 (x,t) + \V 1 (x,t)-2X 2 a 3 ) ^(x,t,X) = 0, 



q(x,t) 
r(x,t) 



CT3 



1 
-1 



where Q(x,t) and 03 are 2r x 2r matrices with compatible block structure. Here 
V 1 (x,t) = 2Q(x,t), V 2 (x,t) = [ad ~}Q,Q] + 2iad ~*Q x (x, t) 



(1.9) 
(1.10) 

(1.11) 
(1.12) 



and ad ~ is the inverse of the adjoint action ad CT3 with respect to the element 03 : ad a3 Y — [(73 , Y] . 

An effective tool to obtain new versions of MNLS is the reduction group introduced by Mikhailov 
|17j . It allows one to impose algebraic constraints on the potential Q(x, t) which are automatically 
compatible with the evolution. For example, the involution (1.6 1, which leads to MNLS with v.b.c. 
(1.7) is known as Z 2 -reduction and can be written as |10j : 

BU\x, t, X*)B~ 1 = U(x, t, X) (1.13) 

where B is an automorphism of g matrix such that B 2 — 1, [03, B] = 0, and 

U(x, t, A) = Q(x, t)-\a 3 . (1.14) 

Below we analyze the multi-component nonlinear Schodinger equation (MNLS): 

wit + q^^Wq + qp + M^o, (lis) 
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with constant boundary conditions (c.b.c.) at x — > ±00: 

lim q(a;,i) = q ± , fjt = q + q\_ = q_q?_ , /I = q+q+ = q^q_ , (1.16) 

where q(x, t) is nxr matrix-valued function, related to A. Ill, C.I or D. Ill-type symmetric spaces. 
The case n r can be related to A. Ill symmetric spaces only and has been solved with the ISM 
in [JT]. 



Therefore we concentrate on the MNLS (1.151 related to C.I or D. Ill-type symmetric spaces, 



which means in particular that n = r. Its Lax pair is obtained from (1.9 1— (1.12) by replacing 
V 2 (x,t) with: 

V 2 (x,t) = [ad-^Q,Q} + 2iad^Q x (x,t)-a 3 Q 2 ± . (1.17) 

Here we have also imposed the additional condition Q\ — Q 2 _ . It ensures that the two asymptotic 
Lax operators L± = + Q± — A03 have the same spectrum. It also ensures that the potentials 
Vi(x,t) and V%(x, t) in the second operator M(A) vanish for x — > ±00. As a result the solutions of 



the MNLS (1.151 q(x, t) do not undergo strong oscillations with respect to time, see [TT1 19"]. 



Lax operators of the form (1.9 1 can be associated with each of the symmetric spaces listed 
below (for the definition see [12] and the Appendix) . They are defined by specifying the simple Lie 
algebra g, having typical representation in 2r x 2r matrices and the Cartan subalgebra element 03: 

• C.I: g ~ C r ~ sp(2r), a 3 — Hg, where the vector a in the root space E r dual to 03 is given 
by 0= ELi e fc 

• D.III : q ~ D r ~ so(2r), cr 3 = H^, where the vector a in the root space E r dual to 173 is 
given by a = J2l=i e k 

Here the orthonormal vectors span the root space E r of both types of algebras. The element 
03 belongs to the Cartan subalgebra f) and is dual to a. Using cr 3 we can split the set of positive 
roots into two two subsets A + = U A±. These sets, for the algebras that we are working with, 
are composed of the following roots: 

A+ = { ei -ej, 1 <i<j<r}, A+ = { 2e 2 , e, + e 3 ■, 1 < i < j < r} (1.18) 

for g ~ sp(2r) and 

A+ = {ei- ej , l<i<j<r}, A+ = { ei + ej , 1 < * < j < r} (1.19) 

for g ~ so(2r). 

The root vectors of the algebra are denoted by E a where a is the corresponding root. 

Let us introduce a projector P a3 = ad~ 3 1 ad CT3 onto the co-adjoint orbit O a3 of the element 03. 
Here the inverse of the adjoint action is ad~ 3 Y = ^u 3 Y. The generic element of A € O ff3 is the 
one that satisfies the relation X = P a3 X. Obviously the potential of the Z-Sh system Q(x, t) and 
its variation SQ(x,t) belong to CT3 . 

This paper extends the results of [TTJ [3] . In Section 2 we focus on the solutions of the direct 
scattering problem for the case of Lax operator describing MNLS with c.b.c. In Section 3 we 
derive the completeness relation for the "squared solutions" of the Lax operator generalizing the 
results of O [TS]. Here we prove that the ISM is equivalent to a generalized Fourier transform 
also for Lax operators with c.b.c. Thus we have shown that the nonlinear evolution of equation 



(1.151 transforms into linear one in terms of the scattering data of L. 
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2 Solutions of the Lax operator L 

The spectrum of the asymptotic operators L± is purely continuous and is determined by the 
eigenvalues of Q± which generically may be arbitrary complex numbers. However, here we consider 
only the case when L becomes self-adjoint. As a result its potential Q(x, t) acquires the form: 



Q(z,i)=-Qt( M ) Q(x,t)=( " A ^ L > (2.1) 



q(x,t) 
-qt(i,t) 

For simplicity reasons we will consider only the case when all of the eigenvalues of the asymptotic 
matrices Q± are real and equal: 

mi = rri2 = ... = m r = m =/= ra£8 (2-2) 

As a result we have the following condition on the eigenvalues of the asymptotic matrices |llj : 
q±q}±(x, t) = m 2 l and the correspondence with the isotropic problem is obvious: fj, = ~p= m 2 l. 

The requirement that the potentials of the Z-Sh system belong to g can be formulated as a 
reduction condition |17l [5] : 

So 1 U t (x,t,\)S = -U(x,t,X), So 1 V t {x,t,\)S = -V(x,t,\), S^^Sq = -a 3 , (2.3) 

which has trivial action on A. The matrix So is the one which realizes the definition of the algebras 
C r ~ sp(2r) or D r ~ so(2r) in the typical representation [5J [T5] . In what follows we will define 
the Lie algebra g by: 

9 = {X:X + S^X'So = 0} , (2.4) 

where 



for g ~ sp{2r) and 



S = ^(-l) s+1 (i? s5 -^ s ) 



5 = ^(-l) s+1 (^ s5 + ^ s ) 



for q ~ so(2r). Here s = 2r — s + 1 and Ek s are 2r x 2r matrices, defined by (Ef~ s )ij = 5ki&sj- Note 
that S 2 = cqI, where eo = — 1 f° r sp(2r) and eo = 1 for so(2r). 

Such reduction (2.3 1 imposes restrictions only on the coefficients of Q(x,t) such that for C r ~ 
sp(2r) we can put: 

r 

Q{x,t) = (qijE ei+ei - q^E_ ei _ ei ) + ^ (q t E 2ei - q*E_ 2ei ) , (2.5) 

i<j i=l 

while in the D r ~ so(2r)-case we have: 

Q(x,t) = J2 (qijE ei+ej - q* t E^ ei ) , (2.6) 

i<j 

where * means complex conjugation. The definitions of the root vectors E a can be found in the 
Appendix. In the typical representations of C r and D r these choices for Q(x,t) have always the 
block structure shown in (2.1|. In the case of g ~ sp(A) the block q is parametrized by three 
functions: _ 
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The corresponding sets of MNLS for these choic es of Q(x, t) and v.b.c. were first derived in [5]. 
For c.b.c. with r = MNLS take the form ( 1.15 1 with the additional linear in q terms ensuring 
regular behavior of the solutions for t — ► ±00. 

Let us outline the construction of the fundamental analytic solutions (FAS). In the particular 
case that we are considering - the isotropic problem - the Jost solutions are defined as fundamental 
solutions with fixed asymptotic for x — » ±00: 

lim xjj(x,X)e ij( - x ^ sX = V>o(A) lim <t>(x, X)e l3{x)a3X = O (A), (2.8) 



x — > — 00 



where 2r x 2r matrices tpo(X) and 4>o(X) take value in the corresponding group Q and diagonalize 
the potential of the Lax operator L 

(Q+ - Aa 3 ) Vo(A) = -^o(A)i(A)a 3 , (Q_ - Aa 3 ) O (A) = -^o(A)j(A)a 3 , (2.9) 
where j(A) = VA 2 — m 2 . They have the block structure: 

*w ^ ( fi i */ ) • *>« - * ( si */ ) • 

The r x r matrices A, B and Si are given by: 



4« = W^^. B H = W^#. S^yiC-l)^.,,-^ (2.11) 



2j(A) ' Ki V 2j(A) 



s=l 



where e P9 are r x r matrices such that (e pq )ij — 5i P Sj q and the phase factor Vo is 2r x 2r diagonal 
and unitary matrix. 

The two Jost solutions arc fundamental solutions and must be linearly dependent. This means 
that there exists a matrix T(t,X), called scattering matrix, which connects them and has an 
appropriate block structure. 

T(t, X) =ip~ 1 (x, t, A) 4>{x, t, A) (2.12) 
The spectral parameter A takes values in two-sheeted Ricmannian surface S: 

S = S 1 L)S 2 

associated with the square root j(X). Each sheet of this surface is determined by the sign of j(A). 

Si : Im j(A) > 0, S 2 : Imj(A) < 0. (2.13) 

Half of the columns of the Jost solutions arc analytic functions of A on the first sheet and the other 
on the second sheet. 

1>(x, A) = (|V>-(z, X)),\iP + (x, A))), cp(x, A) = (\ct>+(x, A)), \<j>~ (x, A))), (2.14) 

where l^^) and |'0 ± ) denote a r x 2r matrix composed of the corresponding r columns of the Jost 
solutions. The superscript " + " means analyticity on the first sheet and "~" - analyticity on the 
second sheet. Next, we can construct FAS on each of the sheets by simply combining the blocks 
of the Jost solutions with the same analyticity properties. 

x +(x,A)= (|0+),|^ + ))(x,A), X -(x,X) = (\^-),\r))(x,X) (2.15) 

Let us write down the FAS x + ( x iX), analytic on the sheet Si and x _ (x, A), analytic on the 
sheet S2 using appropriate decompositions of the scattering matrix (2.191, which consists of the 
same upper (lower) block-triangular functions S 1 * 1 and r T ± as they are in the v.b.c. case [8]: 

X ± (x, A) = ip(x, A)T T = <f>(x, A)S ± (2.16) 
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These triangular factors are: 



and can be viewed also as generalized Gauss decompositions of the T(A). 

T(A) = T~(A)S + (A) = T+(A)S~(A). (2.18) 

Here and after the hat means taking the inverse matrix. We can use for the scattering matrix 
the same block-matrix structure as in v.b.c. case 

«*■ A > " «*■ A)T(A) ' r < A > = ( S(A) "a-w' ) f < A > " ( -d-W ) < 2 ' 19 > 

The elements of the inverse matrix are defined as follows: 

c-(A) = S + (A)(1 + p-p+r 1 = (1 + r+r-)- 1 S + (A) 

d"(A) = a + (A)p"(A) (1 + p+p-y 1 = (1 + r+T-)- 1 r+(A)a-(A) 

c+(A) = a"(A)(l + p+p-y 1 = (1 + r+r-)- 1 a-(A) 

d+(A) - a _ (A)p+(A) (1 + p-p+y 1 = (1 + r-r+)- 1 r-(A)a + (A) 

Here p ± (X) = b ± (A)a ± (A) = B ± (A)d ± (A) and r±(A) = a ± (A)b =F (A) = d =F (A)B ± (A) are the 
multicomponent generalizations of the reflection p ± , t coefficients, (for the scalar case see [2T| 

E2JE9]). 

Given the potential Q(x) one can obtain the Jost solutions uniquely. The Jost solutions in turn 
determine uniquely the scattering matrix T(A) and its inverse T(X). Q{x) contains at most |Ai_| 
independent complex- valued functions of x. Thus it is natural to expect that at most |A^| of the 
coefficients of T(A) for A G R m , instead of (2r) 2 , will be independent. Here |Af | is the number of 
roots in A^,i.e. |A]*~| = r(r+l)/2 for C r and |A^| = r(r — l)/2 for D r . The continuous spectrum 
R m = (— oo, —m) U (m, oo) is determined by the condition |A| > m. 

The set of independent coefficients of T(A) are known as the set of minimal scattering data T. 
They were introduced by Kaup for the Z-Sh system associated with g ~ sl(2) and v.b.c. He proved 
that a (A) can be recovered from T using the analyticity properties, i.e. the so-called dispersion 
relation. The same problem for the generalized Z-Sh system with c.b.c. is more difficult. Here we 
just introduce % — U as the proper generalization of the minimal set of scattering data: 

r 1)C = {p+(A),p-(A), A6U T ltd = {pi(Xf),Xf}f =1 
7i, c = {t+(A), t-(A), A € R m }, r 2>d EE {r±(Xf),Xf}? =1 

where a € A+. The reconstruction of the diagonal blocks a ± (A) from their analyticity properties 
requires a solution of r x r matrix- valued Riemman-Hilbert problem. Here A^ are discrete eigen- 
values of L. The sets % jC characterizing the continuous spectrum need to be completed by the sets 
% y d characterizing the discrete spectrum of L which in turn requires the knowledge of the dressing 
factors. These problems will be addressed elsewhere. 

3 Wronskian relations 

Let the class of allowed potentials M. be a slice of CT3 determined by additional constraints: i.) 
any generic element F(x) = P a3 F{x) of M. is matrix- valued function which vanishes fast enough 



G 



for | a; | — * oo and ii.) the phase factor V which connect the asymptotic values of the potential 
Q + = V^Q-V is an integral of motion. The derivative of the potential Q x (x,t) belongs to the 
class of allowed potentials. The variation of the potential 5Q(x, t) is an allowed potential provided 
it satisfies the second additional condition. The mapping T : M. — > C between the class of allowed 
potentials M. and the scattering data C of L is analyzed by means of Wronskian relations [21 H] • 
These relations allow us to formulate the main result of this work, i.e. that the ISM is a generalized 
Fourier transform in the case of C.I and D. Ill-type symmetric spaces. They also serve to introduce 
the skew-scalar product 

[AOrJ.BOiO] =\ J dx (A(x),[a 3 ,B(x)}) (3.1) 

which is non-degenerate for A(x), B(x) € M. and provides it with symplectic structure. We start 
with the identity: 

(x {Q( x > t) ~ ^3) x(x, A) , E ±a ) I^L.^ = ~i I dx ( -[a 3 ,a 3 Q x ] , P« 3 X E± a x(x, A) 



DC 
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(3.2) 

where x( x ) A) can be any fundamental solution of L. For convenience we choose them to be the 
FAS introduced above. The l.h.side of (3.2) can be calculated explicitly by using the asymptotics 
of FAS for x — > ±00. It would be expressed by the matrix elements of the scattering matrix T(A), 
i.e. by the scattering data of L as follows: 



lP a3 x + (:r, A) E a x + ,a 3 Q x \\ = -j(A) a 3 T~ , E a ) = 2j(A)b+, a e A+ 

lP rT3 X + (x,\)E_ a x + ,<J 3 Q4 =j(\)(s + <r 3 S + ,E_ a )=2j(\)dZ a , aeA+ 

lP <73 X-(x,X)E a x-,a 3 Q4 =j(A)(s~<7 3 S-,£! a )=2j(A)d+ a e A+ 

[P aa x~(x, A) E_ a x~,a 3 QJ = -j(A) (t + ( t 3 T+ , E_ a ) = 2j(X)bZ a , a e A+ (3.3) 



The second set of Wronskian relations which we consider relate the variation of the potential 
5Q to the corresponding variations of the scattering data Sp and St. For this purpose we use the 
identity: 

(x«x(*,A),£?±«)|~_ 00 =y dx ( l -[<T 3 ,a 3 5Q] ,P a3 x(x,X)E ±a xJ (3.4) 

If we assume that the variation of the phase factor 5 V vanishes we arrive at: 
lP a3 X + (x,X)E a x + ,a 3 SQ]] - -i(t~5T-, E, 
lP a3 x + (x,\)E_ a x + ,<j 3 5Q]] =i(s + 5S+,E_ 



= i(5p + a + ) a , 


a G A+ 


-- i((5r+c + )_ Q , 


a G A+ 


= i(6r~c~) a , 


a G A+ 


-- i(6p~a~)- a , 


a G A+ 



\[P a3X -{x,\)E_ aX ;-,a 3 8Q\\=-i{T 5T + ,E_ a \=i(5p-ar)- a , a G A+ (3.5) 

These relations are basic for the analysis of the related NLEE and their Hamiltonian structures. 
The above identities also allow us to introduce the proper generalizations of the usual Fourier 
exponential functions. Let us introduce the set of "squared solutions": 

^t(x,X)=Pa 3 X ± {x,X)E ±a x ± {x,X), foraeA+ (3.6) 
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#£(x, A) = P CT3 X ± (x, X)E Ta x ± (x, A) , for «eA} 



(3.7) 



et{x,X) = P a3 x ± (x,X)E ±aX ± (x,X) 1 foraeA+ (3.8) 
E±{x,\) = P a3 x ± (x,X)E TaX ± (x,X), foraeA+ (3.9) 
r±(x,X) = P rT3 x ± (x,X)H kX ± (x,X) 1 forfc = l,..,r (3.10) 

These are the "squared solutions" of the Lax operator L connected with simple Lie algebra g. 
They are constructed by means of FAS x ± ( x , A) and the Cartan-Weyl basis of the algebra and are 
analytic functions of A on the corresponding sheets of the spectral surface. The equations that 
<I>J an \I>^ satisfy are a direct consequence of the fact that FAS and their inverse satisfy the Z-Sh 
system system: 

+ [Q(x)-\a 3 ,$i(x,\)] = 0, i d ^+[Q(x)-\a 3 ,*t(x,\)]=0 (3.11) 
The " squared solutions" also serve as building blocks of the Green function for L [6l [3 [8] : 

G ± (x, y, A) = Gf(x, y, A)% - x) - Gf(x, y, \)6(x - y), (3.12) 

where 

Gf(x,y,X)= £ $±0r,A)®¥±(y,A) (3.13) 
Gt(x,y,X) = ^2 *£(*,A)®$±(y,A) + £ S± (x, A) ® 0± (y, A) 

a£A| aGA ( |uA + 

r 

+ ^T±(x,A)®T±(y,A) (3.14) 

k=l 

and 0(x) is the usual step function. 



4 Completeness relation and evolution 

The main result in this section is that the sets and {'J J } form complete sets of functions in 

M.. The idea of the proof is simple. Apply the contour integration method along a proper contour 
(see figure|T]) to a conveniently chosen Green function (3.12). From the Cauchy theorem we have 



N 

. f G+(x,y,X)d\ = VRes i G+(x,y,A) (4.1) 

Integrating along the contours we treat separately the contribution from the infinite semi-arcs 
and the ones from the continuous spectrum R m — C\ U C2 which is composed of the cuts C\ = 
(— 00, — m) and C2 — (m, 00). Special care must be taken for the end points A — ±m of the 
spectrum. Assuming that the end points of the spectrum give no contribution 



dX = (4.2) 
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Figure 1: The contour along which we integrate lies completely on the first sheet of the 2-shectcd 
spectral surface associated with the square root j(X) = y/\ 2 — m 2 



where with index e± we have denoted the integrals along the infinitesimal semi-arcs around the 
end points of the spectrum, we obtain the following completeness relation 



5(x - y)U a 



1 T. 

7T ^ 



dX {<f>+(x, A) ® A) - A) ® A)} 



2i EE 



( /A 



a=a: 



A) + A) ® (y, A) 



A=A7 



Here II^ = X^aeA+I^ct ® E- a — -E_ Q ® 25 a ]. The assumption that we have made is that A+ are 
simple poles of the "squared solutions" and 

Using the completeness relation one can expand any generic element of the phase space M. over 
each of the complete sets of "squared solutions" and This relation is utilized with the 
help of the following the trick 



tr x {([as, F{x)} ® l)n CT3 } = itr 2 {n CT3 (l ® [o-a.-Ffar)])} = ^(s) 



(4.3) 



where tr j (and tr 2 ) mean taking the trace of the elements in the first (or in the second) position 
of the tensor product. 

The completeness relation (4.3 1 allows to establish one-to-one correspondence between the 
elements of Ai, such as Q x and Q t , and its expansion coefficients. It is also directly related to 
the spectral decompositions of the generating (recursion) operators A±. These operators are the 
ones whose eigenfunctions are the "squared solutions". Their derivation starts by introducing the 
splitting of the object = x ± (a;, \)E± a x ± (x 7 A) into block diagonal and block off-diagonal parts: 



(x,\) = (l-P a3 )et(x,\) 



(4.4) 
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end use the equation it satisfies 



Thus equation (4.5 1 splits into 



[Q(x),$±(x,A)] =0 



i^+^.ef^A)] =A[(73,*±(a:,A)] 
Equation (4.6 1 can be integrated formally with the result 

e * ± (x, A) = ± (A) + i f dj/ [Q(y) , <i>± (y, A)] 



C*f(A) = lim e^ar.A), 

x — >eoo 



€ = ±1 



Next insert (4.8 1 into (4.7 1 and act on both sides by ad^ 1 . This gives us: 

(A± - A)<&±(x,A) = i [C d a ±{\)M« l 3 Q{x)\ , 



where the generating operators A± are given by 

.dS 



A ± S(x) = ad ff3 {i— 



Q(x) , / dy {Q(y) , E(y)} 



±00 



(4.5) 

(4.6) 
(4.7) 

(4.8) 
(4.9) 

(4.10) 
(4.11) 



Thus (resp. <&J ) will be eigenfunctions of A + (resp. A_) if and only if C^: f (A) = 0. Evaluating 
the limit of (4.9) for all a in the specific case (4.2) we find: 



(A+ - A) A) =0 (A+ - Xf) #±(x, Xf) = 0, 

(A_ - \)$t(x,\) = (A_ - Xf) $t(x,\f) = 0, 

This result can be generalized for arbitrary /(A±): 

(/(A+) - /(A))#±(x,A) = (/(A+) - f(Xf)) *t(x,Xf) = 
(/(A_) - /(A)) $±(x, A) = (/(A_) - f(Xf)) Xf) = 



a € A+ 
a e A+ 



(4.12) 
(4.13) 



a € A+ (4.14) 
a e A+ (4.15) 



The class of higher MNLS on symmetric spaces of C.I and D. Ill-type and with c.b.c. can be 
put down in terms of the derivative of the potential Qt with respect to the evolution parameter 
and the dispersion law /(A) = — 2A [HI IE] as follows: 



iad^Q + Z^ad^Q^O 



(4.16) 



Substituting the objects in this formula with their expansions over the "squared solutions" we 
obtain equations for the evolution of the scattering data. The expansion coefficients of ad^Qt 
and ad ^Q x on the continuous spectrum turn out to be exactly the minimal set of scattering data. 
The evolution for the reflection and transition coefficients is provided by 



.dp ± 



±f(X)j(X)p ± (t,X)=0, 



5r± 



Tf(X)j(X)r ± (t,X)=0 



A e 



(4.17) 
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The observation that the scattering data evolves trivially is visible from the equation depicting 
the evolution of the scattering matrix T(X). This equation is a result of the compatibility condition 
(1.8) and the fact that the t wo Jo st solutions ip and 4> are so lution s of the second operator of the 
Lax pair in the Z-Sh system (1.91. Acting with i4z on T(A) (2.121, we get: 



i~T(t,\)-2\j(\)[a 3 ,T(\)] = 0, 



(4.18) 



where /(A) = — 2A is the dispersion law for the MNLS with c.b.c. This equation for MNLS can 
also be derived from the explicit form of the Lax representation (1.10) by evaluating the limit 
lim-r^-too Mip = 0. For the r x r blocks making up the scattering matrix we have: 



9a±(t,A) 



dh ± 

0, z— T2Aj(A)b ± (i,A)=0 
at 



(4.19) 



These equations have obvious solutions. From the first equation is clear that the diagonal blocks 
are conserved and their invariants - upper (lower) principal minors as well as their determinants 
are generating functionals of the special series of local infinitely many integrals of motion Ik : 



lndeta ± (A) = J2 X ^ Ik 



(4.20) 



k=X 



This is the major idea of the ISM 
q(ai,t), in terms of which the MNLS(1.15l is written 
linear evolution equations. 



a one-to-one change of variables- from the multicomponent 
towards the scattering data which satisfy 



5 Conclusion 

The result of this work is that the interpretation of the ISM as a generalized Fourier transformation 
holds true in the case of Lax operators with constant boundary conditions on symmetric spaces 
connected with the Lie algebras C r ~ sp(2r) and D r ~ so(2r). The completeness relation of 
the "squared solutions" of the generalized Z-Sh system in the case when the Lax operator L 
becomes self-adjoint is derived. The "squared solutions" turn out to be generalizations of the 
usual Fourier exponential function and cigcnfunctions of the recursion operators A±. This result 
allows one to prove that the corresponding NLEE results in linear evolution for the scattering data. 
The recursion operators A± open the path towards the construction of action-angle variables for 
the NLEE solvable with this generalization of the Z-Sh system and from there the Hamiltonian 
formulation of these equations and their hierarchies connected with A±. 

The physical applications of the NLS eq. both with vanishing and non-vanishing boundary 
conditions is well known; the same holds true for the Manakov system as well as for the sp(4) 
MNLS with v.b.c, see [T3]. It will be interesting to find physical applications also for the MNLS 
with c.b.c. 
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Appendix 



The above definition of g (2.4 1 satisfies the requirement that the Cartan subalgebra t) will be made 



up of diagonal matrices. The Cartan generators H^, dual to e^, are given by: 



H, 



E, 



kk 



Ekk 



(5.1) 



The element 03 = X)I=i ^ki belongs to f) and is dual to a. 
The root vectors in the typical representation are given by 



Ej 



E, 



eo(-l) i+J '% 



(5.2) 



where 1 < i < j < r and eo = ±1. Since eo = 1 for g ~ so(2r) equation (5.2) gives vanishing result 
for i = j which is compatible with the fact that 2ei are not roots of so(2r); for g ~ sp(2r) eo = — 1 
and equation (5.2 1 by putting i = j provides also an expression for Ei ev However this expression 

2 



is not normed with respect to the Killing form (E a ,E- c 
with the root 2e^ that we will use are given by [T2"] : 



The Weyl generators associated 



E 9 ^ = V2E, 



(5.3) 
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